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SOME EFFECTS OF COMPRESSIBILITY ON THE FLOW THROUGH FANS AND TURBINES 

By W. Peel and H. T. Epstein 


SUMMARY 

The laws of conservation of mass, momentum, and energy are 
applied to the compressible flow through a two-dimensional 
cascade of airfoils. A fundamental relation between the ulti- 
mate upstream and downstream flow angles, the inlet Mach 
number, and the pressure ratio across the cascade is derived. 
Comparison with the corresponding relation for incompressible 
flow shows large differences. 

The fundamental relation reveals two ranges of flow angles and 
inlet Mach numbers, for which no ideal pressure ratio exists. 
One of these nonideal operating ranges is analogous to a similar 
type in incompressible flow. The other is characteristic only of 
compressible flow. 

The effect of variable axial-flow area is treated. Some impli- 
cations of the basic conservation laws in the case of nonideal 
flow through cascades are discussed. 

INTRODUCTION 


p static pressure 

E gas constant 

S entropy 

s cascade spacing 

t temperature 

u velocity in ^-direction 

v velocity in ^-direction 

w resultant velocity; also mean velocity 

X 2-direction (tangential) force per blade per unit span 

7 ^-direction (normal) force per blade per unit span 

7 ratio of specific beats 

A angle of flow measured from axial direction 
p density 

Subscripts: 

1 upstream of cascade 

2 downstream of cascade 

8 stagnation 


In the endeavor to obtain high pressure rise from axial- 
flow compressors, the operating speeds have been increased 
in recent years to compressibility speeds. The precise extent 
of validity at these speeds of the incompressible-flow methods 
of design of axial-flow compressor blading appears not to 
have been considered heretofore. 

In this paper, which describes results of a theoretical 
investigation made at the NACA Cleveland laboratory during 
1944 and 1945, the extent of validity is determined in the 
isentropic case by applying the compressible-flow laws to the 
ultimate upstream and downstream flow conditions in a 
two-dimensional cascade of airfoils. The resulting equations 
relating the pressure ratio across the cascade, the flow angles, 
and the inlet Mach number are compared with the corres- 
ponding incompressible-flow relations. 

The results apply to both fan-type and turbine-type cas- 
cades. Although compressible-flow relations are generally 
used in the design of turbine cascades and, less generally , 
in the design of fan cascades, it is thought that the compact 
form of the results presented herein will render them more 
easily applicable than the other methods now in use. 

SYMBOLS 

The following symbols are used in this report: 

A flow area perpendicular to axial direction 
c p specific heat at constant pressure 

M Mach number of flow 


VELOCITY-DIAGRAM RELATIONS FOR COMPRESSIBLE 
CASCADE FLOW 

GOVERNING EQUATIONS 

Consider the ultimate upstream and downstream compress- 
ible-flow conditions, or the velocity diagram, for the cascade 
of airfoils representing a fan or a turbine-blade arrangement 
(fig. 1). The conservation laws applied between the up- 
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stream section and the downstream section where the flow 
is assumed to be uniform and to pass through equal areas are: 

Continuity 

Plt![ = p 3 t>2 (1) 

Energy 

( 2 ) 

A A 

Momentum 

7 

Normal, y-direction: — =pi— p 2 +Piffi J — p^ 2 (3) 

s 


For the case of incompressible flow, the expression cor- 
responding to equation (7) may be derived from the equations 
corresponding to (1) and (2); namely 

Continuity (incompressible): Vi=v 2 (8) 

Bernoulli (incompressible): pi+-| PiWi s =JPi+^ piWi 2 (9) 

The result is 

if W 1 -* 1 ®' 1 ) <10 > 


Tangential, ^-direction: ——p\Vi{Ui—vA) ■ (4) 

8 

In addition, there is the gas law 

p= P Rt (5) 

and the isentropic relation 



(The assumption of the isentropic exponent y prevents the 
direct application of the results to fan or turbine blading. It 
appears, however, that use of a suitable polytropic exponent 
n?*y allows direct application in many cases.) 

Substitution of equations (1), (5), and (6) in equation (2) 
yields the fundamental relation between the flow angles Xi and 
X 2 , the pressure ratio PilPi, and the inlet Mach number M\\ 


cos 
cos X 2 




El , 

2 

1 

H 

1 

i 

(■y-1) 


(7) 


in which the substitution 


Hi*- 2 l 

2i 

y~- 

pi 


(ii) 


has been made for comparison with the compressible-flow 
relation (equation (7)). Equation (10) may also be derived 
as the limiting form of equation (7) when the specific-lieat 
ratio y approaches infinity. (The incompressible fluid may 
be regarded as the limiting form of a compressible perfect 


fluid in which the velocity of sound 


V' 


approaches infinity. 


Inasmuch as p and p remain finite, y must become infinite.) 


GRAPHICAL REPRESENTATION OF EQUATIONS 

For a given inlet Mach number M t and pressure ratio 
PtlPu equation (7) shows that the flow angles Xi and X 2 are 
not specified but only their cosine ratio. A plot of the cosine 
ratio cos Xi/cos X 2 against the flow deflection Xj, with 
Xi or X 2 as a parameter (fig. 2), illustrates the possible indi- 
vidual flow angles and deflections for a given cosine ratio. 
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Pressure ratio, pt/p, 

Figure 3. — Relation of compressible flow angles to Pt/Pi for various Inlet Mach numbers. 


The flow angles Xi and Xj may be visualized as the stagger 
angles of the tangents to the camber line of the cascade 
blading at the leading edge and the trailing edge, respectively. 
For low-speed flow, at least, and for shock-free entry, this 
correspondence is approximately valid for cascade solidities 
of the order of unity (solidity=blade chord/blade spacing) 
and increases in validity with the solidity. 

The contours representing equation (7) are shown in figures 
3 and 4. Several of the contours of figures 3 and 4 are com- 
pared with the corresponding contours in the incompressible 
case (equation (10)) in figures 4 and 5. 

COMPARISON OF COMPRESSIBLE AND INCOMPRESSIBLE EQUATIONS 

The comparison between the compressible and incompres- 
sible equations (figs. 4 and 5) reveals considerable differences 
except in the pressure-rise range (p 2 /?i))>l at low inlet Mach 
number (Af : <0.4). Outside this range consider, for exam- 
ple, the determination of the blade camber of a high-solidity 
pressure-rise cascade designed to operate at an inlet-flow 
angle of 45°, an inlet Mack number of 0.80, and a pressure 
ratio of 1.2. (Camber is defined as the difference in angle 
between leading-edge and trailing-edge camber-line tangents. 
Camber therefore approximates X t — Xj for high-solidity 
cascades.) By the incompressible contour 11/1=0.8 in figure 
5 and by use of figure 2, the camber should be 26.9°. Accord- 
ing to the compressible curve, on the other hand, the camber 
should be 9.4°, or 65 percent smaller. If the inlet Mach 
number is increased to unity, the results are: 


Incompressible camber, degrees 1L 8 

Compressible camber, degrees L 5 


As an example of the differences between compressible and 


incompressible flow in the case of a pressure drop through the 
cascade, assume X!=0°, A/ =0.80, and p*(pi=0.7. Then 


from figures 2 and 5 

Incompressible-flow deflection, degrees .. 39. 3 

Compressible-flow deflection, degrees 13. 1 


Thus, the effect of compressibility in the normal operating 
range of moderate pressure ratios is to lower substantially 
the required flow deflection for a given pressure ratio. 

On the other hand, for given flow angles, that is, given 
cos Xi/cos Xj, figure 4 shows that for normal compressor 
operation, (cos X 1 /cosX 2 )<Cl and (ps/^i)>l, the compressible 
pressure ratio increases with the inlet Mach number Mi 
at a greater rate than the incompressible pressure ratio. 
For example, if cos Xi/cos X2=0.8 in normal compressor 
operation and ilfi=0.8, the incompressible pressure ratio is 
1.16 and the compressible pressure ratio is 1.27. 

FURTHER DISCUSSION OF THE COMPRESSIBLE CASE 

The surface representing the fundamental equation (7) is 
the locus of all possible isentropic {loss-free) operating 
conditions of the cascade and no others. The contour projec- 
tions of this surface show that certain ranges of operating^ 
conditions are unattainable in isentropic flow; thus, in figure 
5, there is no inlet Mach number Mi to correspond to cosine 
ratios cos Xj cos Xj and pressure ratios ptjpi in the plane 
regions AOB and COD. 

The reason for such nonideal regions can be seen by analogy 
with the incompressible case. The incompressible M\ con- 
tours entirely exclude the first and the third quadrants, 
relative to point O as origin. Consider the case 



. .. COS Ai 

Cosine ratio, CO g j Mach number 
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Figure 5. — Comparison of compressible and incompressible flow angles as function of psjpi for various inlet Mach numbers. 
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(cos Xi/cos Aj}> 1 . Inasmuch as the exit angle Xa is greater than 
the entrance angle X x , the exit tangential-velocity component 
is greater than the entrance tangential-velocity component. 
Because the axial velocity is constant, a pressure drop 
(Ps/Pi)^! is required in accordance with Bernoulli's 
theorem, and the operation point for (cos Xi/eos X0>1 cannot 
lie in the first quadrant but must lie in the second quadrant . 
Similarly, operating points for (cos Xi/cos \ 2 )>1 must lie 
in the fourth quadrant. The partial overlapping of the 
compressible-flow curves into the first and third quadrants 
is due to the influence of compressibility on the axial com- 
ponent of velocity. 

Another type of nonideal region is shown in figure 4. For 
fixed (cos Xi/cos X 2 )>1, an isentropically impossible range 
of inlet Mach numbers exists, which increases with cos Xi/cos X s . 
At cos Xi/cos Xa = 1.1, for example, the excluded 
Mi range is 0.695<Mi<1.37. This type of nonideal range 
is not present in the incompressible case shown as the 
dashed curves of figure 4. The lower limiting inlet Mach 
numbers of this range correspond to maximum flow through 
the cascade, for, because cos Xi/cos X 2 is constant along each 
curve in figure 4, the maxima of M x also correspond to the 
maxima of M x (cos Xi/cos X s ). For a given high-solidity 
cascade, Xj varies only slightly as the inlet conditions are 
varied; therefore, M x cos X 1( which is proportional to the 
flow through the cascade, is likewise a maximum at these 
points. This particular limiting Mach number is indepen- 
dent of the blade-shape details and, hence, of the ratio of 
minimum flow area to inlet area. 

If a cascade is designed to operate in an isentropically 
impossible region, presumably much greater losses, caused 
by flow separation and shock waves, will occur than would 
be encountered at an ideally permissible operating point. 
Such nonideal designs appear to be not at all unusual on the 
basis of incompressible-flow theory. 

Of considerable interest is the fact that, at the maxima 
and the minima of all cos Xi/cos Xa contours, the exit Mach 
number M 2 is unity, which may be shown as follows: The 
maxima and the minima of the cos Xi/cos X 2 contours in figure 
4 correspond to the maxima of the M x contours in figure 3. 
The loci of these maxima are indicated in both figures. The 
condition for maximum cos Xi/cos X* along an M x contour is 
obtained by differentiation of equation (7) as 

Af, cos Xl =(2lW ( 12 ) 

1 cos X 2 \j?i/ 


The downstream Mach number 


A/,= 


W£_ 

Pi 
J — 
Pi 


(13) 


Substitution of equation (12) in equation (14) yields the 
desired result, namely, M*=l. In the limiting case Mi=0, 
equations (12) and (7) show that 

/ cos XA = ^ 

\eos X 2 /M 1-0 



(y=1.4) 


This limiting pressure ratio is the same as the critical 
pressure ratio for local Mach number of unity in one- 
dimensional compressible flow. 

Points to the right of the maximum on an M r contour 
(fig. 3) correspond to subsonic exit conditions M*<1; 
whereas points to the left of the maximum correspond to 
supersonic exit conditions M a > 1. In figure 4 the subsonic 
exit region is to the right of the locus of maxima and minima 
and the supersonic exit region to the left. The proof is as 
follows : From general principles, the ideal continuous flow 
through a cascade is reversible; that is, equation (7) remains 
valid if the subscripts 1 and 2 are interchanged. This 
validity may also be shown directly by substitution of 
equation (14) in equation (7) to eliminate M x . Consequently, 
the exit Mach number M a can be determined from the 
curves that represent equation (7); thus, in figure 3, the exit 
Mach number M 2 for the operating conditions M u p 2 /p X) and 
cos Xi/cos X 2 is obtained by interpolation from the M contours 
at the point for which the coordinates are the reciprocal 
quantities pilp 2 and cos Xj/cos \ x . This reciprocal relation, 
applied to the loci of the maximum points on the M t contours 
(fig. 3), yields the M x —1 contour (where now M x =M t ), by 
the proof previously given; therefore, the reciprocal relation 
applied to a point to the right of a maximum yields a point 
in the subsonic region relative to the M x =l con torn 1 and 
vice versa. 


VARIABLE AXIAL-FLOW AREA 

In many multistage compressors and turbines, the axial- 
flow area is continuously varied in order to maintain the 
axial-flow velocity constant. The effect of a change in 
axial-flow area through a cascade can be very simply taken 
into account in the preceding analysis. The equation of 
continuity (equation (1)) is now 


or 


p x v x A x —piViAi 

(15) 

p x v x a=p 2 v 2 

(16) 


where A x and A 2 are the axial-flow areas upstream and down- 
stream, respectively, of the cascade, and 


can now be written, with the aid of equations (1), (2), (5), 
and (6), as 


M 2 


M cosXi 
cos X 3 



(14) 



(17) 


The only change required in the fundamental equation (7) 
is that the flow-angle parameter cos \ x f cos X s be replaced by 
a (cos X^cos X 2 ). 
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ENERGY RELATIONS FOR CASCADES 

If the isen tropic relation (equation (6)) is dropped, equa- 
tions (1) to (5) relate ultimate conditions on each side of a 
cascade within which losses may occur. The equations may 
also be regarded as relating conditions across two sides of a 
line of discontinuity separating two regions of uniform flow, 
with the force densities Xjs, Y/s acting on the fluid along the 
discontinuity. The oblique shock wave in supersonic flow 
represents the special case in which X/s=Y/s=0. 

Without the isentropic relation there is, of course, no longer 
a unique dependence of the flow-angle ratio cos \Jcos X 2 
on the pressure ratio p 2 /pi and the inlet Mach number Mi. 
(See equation (7).) Equations (1) to (5), however, yield 
information concerning the relation between the pressure 
ratios that might actually occur, the forces on the blade, 
and the losses in the cascade. Suppose, for example, that 
a high-solidity compressor cascade were designed for a certain 
ideal operating point. ' The ideal forces on the blades could 
be calculated from equations (3) and (4). In the actual 
flow the pressure ratio would be less than the ideal value. 
The flow angles, however, because of the high solidity, could 
be expected to remain nearly unaltered. On the assumption 
of an ideal flow-angle ratio and a smaller than ideal pressure 
ratio, the actual blade forces could be computed by equa- 
tions (3) and (4). The vector difference between the actual 
resultant force and the ideal resultant force is a difference 
force caused by friction and shock-wave losses. The corre- 
sponding stagnation-pressure loss or, equivalently, the 
entropy increase across the cascade, can also be calculated. 

An example of the type of computation just outlined fol- 
lows: Among other tilings, it will be shown that the difference 
force is neither at right angles to the ideal resultant force nor 
in the direction of the mean relative flow. Although, as 
will be shown, the incompressible equations predict that the 
difference force is always in the axial direction, the compressi- 
ble equations show that this is far from being the case. 

The losses in cascade flow between an ultimate upstream 
and an ultimate downstream uniform-flow condition may be 
expressed in terms of the actual static-pressure ratio pt/pi 


and the density ratio p a /pi as follows: In the downstream 
condition, the ratio of stagnation pressure p,,i to static pres- 
sure Pi is isentropically related to the corresponding tempera- 
ture ratio as 

/ Pt.iYY 
\PiJ U 

Similarly, in the upstream condition, 

(Vi^VTJ-jA 
\pj k 

If the preceding two equations are^divided, remembering 
that because no energy is supplied from the outside, 

and if pressure and density are substituted for temperature 
by the gas law, there results 


Pa. 


(pi 

APi 


y 

\r— l 


J 


PVL 


(18) 




In isentropic flow, equation (18) reduces, of course, to 
p,.i=p,.\. In an actual flow p li2 is less than p t .\ corres- 
ponding to an entropy increase Si-Si, given ;by the well- 
known formula 


Si-S^R log. ^ 


(19) 


Substituting the normal velocity ratio r 2 /»i into equation (18) 
yields, by the continuity equation (1), 



The normal velocity ratio can bo expressed in terms of 
the static-pressure ratio p%lp\, the inlet Mach number Af lt 
and the cosine ratio cos X t /cos X 2 by equations (1), (2), and 
(5). The result is 


r s 


i Pi 


LV 


1 + 2 (y— 1) Mi 2 

w COS 2 X 2 


(7—1) 


(fj 




( 21 ) 


also 


2 Vi t’i tan X 2 
u~ Vi uC~ tanXi 


( 22 ) 


The expressions (3) and (4) for the blade forces may be 
written in nondimensional form as 


-L =1 _g_ 7Ml wx i 


GrO 


(23) 


§r yM * cosS Xi 0 


tan X a 
tan Xi vj 


(24) 


As examples of the use of the preceding equations, losses 
aDd blade forces for a pressure-rise and a pressure-drop cas- 
cade have been computed and the results are listed in table I. 
For both cases an inlet Mach number of unity was assumed 
with inflow and outflow angles of 55° and 45°, respectively. 
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These flow angles, though possible for pressure rise, are im- 
possible for pressure drop by incompressible-flow theory 
(under the condition of constant axial-flow area assumed 
here). The ideal quantities were first calculated and then 
for the loss case, new pressure ratios were assumed, which 
yielded the stagnation-pressure losses shown. It will be 
noted that the static-pressure ratio of 0.5, which yields a 
stagnation -pressure ratio of 0.942 in the pressure-drop case, 
is greater than the ideal static-pressure ratio. In incom- 
pressible flow it would have to be lower. 



Fionas 6.— Typical cascade forces and velocities. 

The cascade-flow angles and forces are illustrated in figure 
6. The resultant force l, the difference force d, and the mean 
velocity w bear no simple directional relation to each other. 
Because, in incompressible flow (under the conditions of 


constant axial-flow area assumed here), t>a/i>i=l and, in high- 
solidity cascades, Utju-i remains unchanged from the ideal case 
by equation (22), it is evident from equations (3) and (4) 
that only the static pressure can change in going from the 
ideal to the loss case. The incompressible difference force 
must therefore be in the axial direction. This conclusion 
for the incompressible case was reached in reference 1 by 
another method. The losses in high-speed cascade flow must 
evidently be considerably influenced by compressibility 
effects. 

CONCLUSIONS 

1. For inlet Mach numbers less than 0.4 in the normal 
operating range of an axial-flow pressure-rise stage, the 
relations for compressible and incompressible flow yield 
very nearly the same results. For larger inlet Mach numbers 
the discrepancies become considerable. 

2. The compressible equations show the existence of 
nonideal operating ranges, which are not indicated by the 
incompressible equations. Conversely, certain operating 
ranges, which are excluded by the incompressible equations, 
are possible according to the compressible equations. 

3. In certain ranges of the flow angles, two possible pres- 
sure ratios are predicted for given inlet Mach numbers and 
flow angles. 


Aircraft Engine Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, August 1, 1945. 
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TABLE I- LOSSES AND BLADE FORCES FOR CASCADES IN 
COMPRESSIBLE FLOW 


[X 1 =55°;X I =45°;M 1 =1] 



Pressure rise 

Pressure drop 

Ideal 

With 

losses 

Ideal 

With 

losses 

Pressure ratio, pi/pi — 

Stagnation pressure ratio, 

Tangential force, X/t pi 

Normal force, 37* Pi 

Eatlo of resultant force to difference 

1.522 

1 

0-315 
-a 403 

00 

*1.4 

aflso 

0- 290 
-a 307 

5.17 

0.«1 

1 

-0.113 

0.19# 

a* 

*0-5 

0.942 

- 0.111 

0.192 

7.60 



» Assumed values. 




